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CHAPTER  I 
INTRODUCTION 


This  chapter  contains  a  statement  of  the  problem  considered  in  the 
dissertation,  the  motivation  for  considering  the  problem,  a  summary  of 
the  contents,  and  a  brief  survey  of  the  literature  on  the  control  of 
electromagnetic  scattering  from  radar  targets. 

The  objective  of  the  research  reported  herein  is  to  investigate 
the  application  of  impedance  loading  techniques  for  the  reduction  of 
backscatter  from  radar  targets.  A  primary  goal  is  to  achieve  tech¬ 
niques  for  minimizing  backscatter  by  controlling  edge  diffraction 
sources.  Edges  are  major  contributors  to  backscattering  for  many 
scattering  geometries  such  as  wing  structures  and  finite  conical 
shapes. 

The  suppression  of  scattering  is  a  camouflage  problem,  but  over 
the  years  the  development  of  techniques  for  the  reduction  of  specular 
contributions  has  received  most  attention.  Although  the  specular 
contributions  tend  to  be  restricted  in  the  aspect  angles  for  which  they 
occur,  they  dominate  the  net  return  whenever  they  are  present.  How¬ 
ever,  the  specular  echoes  are  amenable  to  suppression  using  radar  ab¬ 
sorbers.  On  the  other  hand,  non-specular  scattering  is  more  pervasive. 
Even  if  all  specular  echoes  are  suppressed,  the  non-specular  scattering 
still  can  generate  radar  returns  of  unacceptable  magnitude.  In  fact, 
modern  aerospace  vehicles  are  often  not  seen  at  specular  aspects  in  a 
tactical  environment.  But  the  non-specular  radar  cross  section  can 
still  reach  levels  that  increase  the  probability  of  detection.  It  is 
therefore  important  to  study  the  backscatter  reduction  for  non-specular 
aspects. 

In  Chapter  II  the  impedance  boundary  condition  and  the  integral 
equations  appropriate  in  an  antenna  and  scattering  problem  are  ex¬ 
amined.  The  low  frequency  numerical  technique,  the  method  of  moments, 
and  the  high  frequency  wedge  diffraction  theory  are  also  discussed. 
These  theories  provide  the  necessary  theoretical  background  for  the 
applications  in  the  later  chapters. 

Chapter  III  presents  the  hybrid  methods  which  combine  the 
Geometrical  Theory  of  Diffraction  (GTD)  and  the  Method  of  Moments. 

These  hybrid  methods  can  be  used  to  treat  problems  which  can  not  be 
solved  by  either  GTD  or  moment  method  alone. 
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In  Chapter  IV,  the  hybrid  methods  are  applied  to  study  the  back- 
scatter  reduction  from  various  geometries  such  as  the  infinite  wedge, 
the  thin  strip  and  a  two-dimensional  wing  model.  The  loadings  required 
for  optimum  reduction  are  also  shown. 

In  Chapter  V,  the  scattering  from  a  thin  square  plate  loaded  with 
a  slot  has  been  studied  as  a  way  of  implementing  the  required  antenna 
impedance  loading.  The  effect  of  loading  as  a  function  of  frequency  is 
also  presented. 

A  considerable  amount  of  effort  has  been  devoted  to  the  reduction 
of  backscatter  from  various  radar  target  geometries  [1 ,2,3,4]  over  the 
years.  For  modifying  the  radar  cross  section  of  scattering  bodies, 
two  techniques  have  been  used  either  separately  or  jointly;  shaping  of 
the  body  and  coating  of  the  body  with  radar  absorber.  Recently,  the 
surface  impedance  model  of  absorber  and  the  integral  equations  were 
used  by  Knott  et  al .  [5,6,7]  to  study  the  nor-specular  radar  cross 
section. 

A  more  encouraging  method  of  reducing  radar  cross  section  is  the 
antenna  impedance  loading  techniques  [8,9].  The  scattering  from  a 
conducting  object  can  be  controlled  by  loading  portions  of  its 
surface  with  distributed  or  lumped  impedances.  The  scattered  field 
from  a  loaded  object  can  be  expressed  as  a  superposition  of  fields 
scattered  by  the  unloaded  object  and  fields  reradiated  by  the  loading. 

The  interference  (  these  fields  can  therefore  be  used  to  control  the 
scattering.  The  impedance  loading  technique  has  been  applied  to  dipole, 
cylindrical  and  loop  antennas  to  minimize  the  radar  cross  section  [10,31]. 
The  backscatter  reduction  in  the  broadside  direction  for  a  small  plate 
loaded  with  an  open  or  short  circuited  slot  was  investigated  by 
Green  [32]. 


CHAPTER  II 

THEORETICAL  BACKGROUND 


A.  Introduction 

In  this  chapter  the  impedance  boundary  condition  and  the  integral 
equations  appropriate  in  an  antenna  and  scattering  problem  will  be 
examined.  The  purpose  of  this  study  is  to  show  how  the  equations 
are  affected  by  the  geometry  and  the  boundary  conditions.  A  general 
procedure  for  solving  linear  equations,  the  method  of  moments,  will 
be  discussed.  The  wedge  diffraction  solution,  which  is  a  very  im¬ 
portant  part  of  the  Geometrical  Theory  of  Diffraction,  is  also 
presented.  Each  of  these  solutions  are  presented  in  their  basic 
terms  in  this  chapter  and  applied  to  specific  problems  in  the  fol¬ 
lowing  chapters. 


The  Impedance  Boundary  Conditions 


In  its  most  straightforward  form  an  impedance  boundary  con¬ 
dition  [11]  is  one  which  relates  the  tangential  components  of  the 
electric  and  magnetic  fields  via  an  impedance  factor.  The  impedance 
factor  itself  is  a  function  of  the  properties  of  the  surface.  The  con¬ 
ditions  are  widely  used  in  diffraction  problems  and  the  analysis  of 
surface  wave  phenomena.  In  this  dissertation,  the  surface  impedance  is 
employed  to  represent  a  loaded  slot.  Consider  a  closed  cylindrical 
surface  whose  generators  are  parallel  to  the  z  axis  of  a  cylindrical 
coordinate  system  and  whose  profile  in  a  plane  perpendicular  to  the  z 
axis  is  the  closed  curve  c.  Let  n  be  an  outward  unit  vector  to  c  and  let 
t  be  a  unit  vector  in  the  tangential  direction,  such  that  t,  n  and  z  form 
a  right-hand  system.  Then,  on  the  surface,  the  impedance  boundary  con¬ 
dition  is  given  by  [12,13] 

■  ■  .  A  A  _  A 

(1)  E  -  (n-E)n  =  Z$  n  x  H 


where  E  and  H  are  the  total  electric  and  magnetic  field,  respectively, 
in  the  region  surrounding  the  body,  which  is  that  of  free  space  for 
the  radar  environment  and  Zs  is  the  surface  impedance.  Zs  may  be  a 
function  of  distance  t  along  the  cylindrical  surface  c.  For  a  perfectly 
conducting  surface,  Zs  =  0.  Equation  (1)  can  alternatively  be  ex¬ 
pressed  as  [12,13] 

V  .  A  -  V  A  1  A  «. X 

(2)  H  -  (n-H)n  =  -  j-  n  x  E 

Ls 

Equation  (2)  is  the  dual  of  Eq.  (1)  and  is  not  valid  if  at  any  point 
on  the  surface  Z$  is  zero. 
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C.  The  Reaction  .ntegral  Equation 


The  reaction  concept  and  its  applications  have  been  discussed  by 
Rumsey  [14],  Cohen  [15],  Harrington  [16]  and  Richmond  [17].  The 
reaction  integral  equation  has  been  used  extensively  by  Richmond  [19, 
30,33]. 

Let  us  consider  the  application  of  Schelkunoff 's  surface  equiva¬ 
lence  principle  [18]  to  the  exterior  scattering  problem  illustrated  in 
Fig.  la.  In  the  presence  of  a  dielectricjorjsonducting  uody,  the 
external  electric  and  magnetic  c^rr§nts  (J-j  ,M-j )  generate  the  electric 
and  magnetic  field  intensities  (E,  H).  For  convenience,  assume  the 
exterior  medium  is  free  space. 

According  to  the  equivalence  principle,  if  we  introduce  the  fol¬ 
lowing  surface  current  densities 

(3)  =  n  x  H 

(4)  =  E  x  n 


on  the  closed  surface  S  cf  the  scatterer,  the  interior  field  will 
vanish,  without  disturbing  the  exterior  field.  Since  the  field  in¬ 
ternal  to  S  is  zero,  as  illustrated  in  Fig.  lb,  we  can  replace  the 
scatterer  by  free  space  without  cjist^urbing  the  field  anywhere.  By 
definiti^n^the  incident  field  (E^.H-j)  is  generated  by  the  external 
source  in  free  space,  and  the  scattered  field  is: 


(6)  H s  =  H  -  H i  . 


From  the  superposition  prin,cigle  and  the  uniqueness  theorem,  we  know 
that  the  surface  current,  (JS,MS)  which  radiates  in^fre^  space  will 
generate  the  f  ield  (ES,HS)  in  the  exterior  and  (-E..,-^. )  in  the 
interior  region,  as  shown  in  Fig.  lc. 
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(E,H) 


Fig.  la— Tj}e  ^source  (J^  ,M-j )  generates  the  field 
(E,H)  with  scatterer. 


Fig.  lb--The  interior  field  vanishes  when  the  current  (JS,MS) 
are  introduced  on  the  surface  of  the  scatterer. 


,C>s 

\  N 

\  \ 


V\ 


FREE  SPACER 


Fig.  1  c— Tlie  exterior  scattered  field  may  be  generated  by 
(JS,MS)  in  free  space. 
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With  the  scatterer  replaced  by  free  space  and  the  equivalent 
current  on  its  surface  S,  we  have  noted  in  Fig.  lb  that  the  interior 
region  has  a  zero  field.  Now  let  us  place  a  test  source  or  probe  in 
the  interior  region,  as  shown  in  Fig.  2,  and  consider  its  reaction  with 
the  other  sources.  If  the  test  source  has  electric  current  density  J*. 
and  magnetic  current  density-^,  the  reaction  principle  gives 


(7) 


II  <ff 


Es-Mt.Hs)ds  = 


<aVEVvH> 


where  (ES,HS)  is  the  scattered  field  generated  by  (Cfs,tfs),  and  the 
integrals  extend  over  the  surface  of  the  test  source.  Equation  (7) 
is  one  form  of  the  reaction  integral  equation  (RIE).  As  noted  by 
Richmond  [19],  the  RIE  is  more  general  than  the  electric  field  integral 
equation  (EFIE)  or  the  magnetic  field  integral  equation  (MFIE).  If  we 
enforce  Eq.  (7)  with  a  set  of  delta -function  electric  test  sources,  the 
RIE  reduces  to  the  EFIE.  If  we  enforce  Eq.  (7)  with  a  set  of  delta- 
function  magnetic  test  sources,  the  RIE  reduces  to  the  MFIE. 

From  Eq.  (7)  and  the  reciprocity  theorem,  we  obtain  another  form 
of  the  reaction  integral  equation: 


(8) 


<fj>  (tj'E't  -  pfs'fit)ds  +  /If  (;Vft  ‘  "i,Fit)dv  =  0 


where  (tt»lft)  is  the  field  of  the  test  source  radiat.ng  in  free-space. 
In  words,  Eq.  (8)  states  that  the  interior  test  source  has  zero  re¬ 
action  with  the  other  sources.  This  zero-reaction  was  developed  by 
Rumsey  [14], 

Combining  the  impedance  boundary  condition,  Eq.  (1),  with  the 
equivalent  current  equations,  Eqs.  (3)  and  (4),  gives: 


(9)  FTS  =  Zs  X  S  . 

Equation  (9)  is  another  form  of  the  impedance  boundary  condition.  For 
two-dimensional  problems  involving  cylindrical  scatterers,  Js  and 
are  functions  only  of  the  position  t  around  the  cylindrical  surface  c. 
If  J.  vanishes,  Eqs.  (8)  and  (9)  yield 

(10)  •  [£t  -  (n  x  Ht)  Z$]dt  =  jj  •  tft  ds. 
c 
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Fig.  2  A  test  source  is  placed  in  the  interior  region 

of  the  scatterer. 


D.  Method  of  Moments 

The  method  of  moments  [20]  is  a  general  procedure  for  solving 
integral  equations.  It  is  widely  used,  with  the  help  of  high  speed 
digital  computers,  in  solving  low-frequency  electromagnetic  antenna 
and  scattering  problems.  Consider  the  linear  inhomogeneous  equation 

01)  Lop(J).E 

where  Lop  is^a  linear  operator,  E  is  a  known  excitation  function  or 

source,  and  J  is  the  unknown  response  function  to  be  determined. 
Equation  (11)  may  be  the  RIE,  EFIE  or  MFIE. 

The  procedure  for  obtaining  an  approximate  solution  to  Eq.  (11) 
can  be  summarized  as  below.  First,  expand  the  unknown  function  in 
a  series  of  basis  functions  in  the  domain  of  the  operator  LQp,  as 

(12)  ^ 
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Substituting  Eq.  (12)  into  (11)  gi>es 


The  next  step  is  to  define  a  set  of  weighting  functions  w, ,  •••  in 

the  domain  of  LQp  and  then  form  the  inner  product  L 


(14) 


a  <  w.  L  J  >  =  <  w  .  E>. 
n  m  op  n  m 


Equation  (14)  is  a  set  of  linear  equations  which  can  be  solved  for  the 
unknown  ap  by  numerical  techniques. 

Two  special  cases  of  the  moment  method  are  of  particular  im¬ 
portance  in  electromagnetics,  namely;  the  Galerkin's  method  and  the 
point  matching  technique.  The  Galerkin's  method  uses  the  same  set  of 
functions  as  the  expansion  and  weighting  functions,  It  requires  the 
evaluation  of  integrals  which  may  be  tedious  and  computationally  ex¬ 
pensive.  This  difficulty  can  be  minimized  if  one  uses  point  matching 
which  uses  the  Dirac  delta  functions  as  testing  functions. 

E.  Wedge  Diffraction 

Wedge  diffraction  is  a  useful  tool  in  the  analyses  of  high- 
frequency  electromagnetic  problems.  The  asymptotic  solution  for 
the  diffraction  from  a  conduting  wedge  was  first  obtained  by 
Somme rfeld  [21].  Originally,  GTD  [34]  as  applied  to  diffraction  by  a 
wedge  was  based  on  plane  wave  diffraction  coefficients;  however,  as 
shown  by  Russo,  Rudduck  and  Peters  in  Refs.  [22,41],  the  use  of  dif¬ 
fraction  of  cylindrical  waves  has  been  found  necessary  in  the  treatment 
of  antennas.  Consequently,  different  formulations  of  wedge  diffraction 
were  substituted  for  the  plane  wave  diffraction  coefficient  which  is  the 
basis  for  wedge  diffraction  theory.  Pauli  [23]  introduced  the  Vg 
function  as  ?  practical  formulation  to  the  solution  for  a  finite-anglp 
conducting  wedge.  Recently,  however,  Hutchins  and  Kouyoumjian[24, 

25]  presented  a  formulation  for  the  diffracted  field  (Vg),  which  sig¬ 
nificantly  improves  the  accuracy  over  that  obtained  from  Pauli's 
form. 


This  improved  diffraction  solution  [24,25]  is  better  in  the 
transition  regions  (near  the  incident  and  refleted  shadow  boundaries). 
It  can  be  written  in  the  form 

(15)  VB(L,e,n)  =  I^L.e.n)  +  I+7T(L,e,n) 
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where 


-j(kL+„/4) 

I  (L,3,n)~  - — - Ja  cot 

jni27 


TT±  3 

2n 


x 


ejkLa 


dx  +  [higher  order  terms] 


where  the  higher  order  terms  are  negligible  for  large  kr  and 
with  n  defined  from  the  wedge  angle  WA  =(2-n)n,  also  a  =  1  +  cos(e  - 
2nirN)  and  N  is  a  positive  or  negative  integer  or  zero,  whichever  most 
nearly  satisfies  the  equations 


2n-rrN  -  e  =  -it  for  I 


2mrN  -  3  =  +tt  for  I, 

I'll 

Consider  a  plane  electromagnetic  wave  normally  incident  on  a  wedge 
of  angle  (2-n)Tr  as  illustrated  in  Fig.  3.  Cylindrical  coordinates 
are  employed  here  with  the  z-axis  coincident  with  the  edge  and  con¬ 
sequently  normal  to  the  plane  of  diffraction.  The  z-component  of  the 
total  field  at  observation  point  P(r,ijj)  may  be  represented  by  the 
scalar  function  u,  given  by 

(16)  u(r,ifi,n)  =  V(r,(M,0,n)  -  V(r,^o,n)  . 

The  minus  sign  applies  for  the  boundary  condition  u  =  0  on  the  wedge 
surface;  thus  the  solution  is  valid  for  the  electric  field  oriented 
in  the  z-direction  (parallel  to  the  edge).  The  plus  sign  applies 
for  the  normal  derivative  of  u  equal  zero  on  the  wedge  surface,  or 
the  magnetic  field  oriented  in  the  z-direction  (polarization  per¬ 
pendicular  to  the  edge).  The  quantities  V(r,4>,n)  represent  the 
separation  of  the  solution  into  incident  and  reflected  waves  and 
are  given  by 


(17)  V(r,<t>,n)  =  V‘(r,(j>,n)  +  Vgfr^.n) 

where  V*  is  the  geometrical  optics  wave  which  is  given  by 


9 


8 


V*(r,4>,n) 


r 


exp[jkr  cos(<j>+2irnN)],  if  -Tr«t>+2TrnN<TT 

for  N=0,  ±1  ,±2  •  •  • 


0,  otherwise. 


Fig.  3--Diffraction  of  a  plane  wave  by  a  conducting  wedge. 
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CHAPTER  III 
HYBRID  METHODS 


A.  introduction 


For  many  years,  GTD  has  been  applied  to  high-frequency  antenna 
and  scattering  problems  for  which  the  structure  is  large  in  terms  of 
wavelength.  GTD  solutions  have  mainly  been  used  in  pattern  com¬ 
putations  with  little  information  concerning  antenna  impedance 
and/or  aperture  distribution.  This  is  especially  true  for  wire- 
type  antennas.  On  the  other  hand,  the  method  of  moment  has  been 
applied  to  low-frequency  problems  for  which  the  structure  is  small 
in  terms  of  wavelength. 

A  hybrid  method  (MMGTD)  to  combine  the  GTD  and  thf'  Moment 
Method  was  used  by  Burnside  [27,28]  to  investigate  the  scattering 
of  a  TE  plane  wave  by  infinite  conducting  structures.  In  this 
chapter,  the  original  hybrid  method  is  extended  to  treat  infinite 
geometries  with  local  inhomogenities  which  can  not  be  solved  by 
either  GTD  or  moment  method  alone.  While  the  method  gives  good 
results  for  finite  geometries,  it  fails  to  predict  the  fields  of 
infinite  structures  accurately.  A  modified  hybrid  method,  which  is 
more  accurate  and  economical  in  computation,  is  introduced.  The  modi¬ 
fied  method  can  be  applied  to  geometries  with  local  surface  impedance 
if  the  surface  current  density  for  the  conducting  case  is  known.  The 
modified  method  was  first  suggested  by  Rudduck. 


B.  A  Hybrid  Method 


Consider  the  scattering  by  a  simple  perfectly  conducting  wedge  for 
a  TE  plane  wave  incident  as  shown  in  Fig.  4.  At  the  field  point  (p,i i>), 
the  incident  field  is  given  by 


jkP  cos(^-^0) 
e 


for  |  <  18°°  . 


The  reflected  field  from  the  wall  along  the  x-axis  of  Fig.  4  is  given 
by 


jkp  cost  ft* ) 
lf2(p.*)  =  e  0 


for  j^+i^0|  <  180° 


The  diffracted  field  can  be  expressed  as  [26] 
H^(p>^)  =  Vg(p,«p-\p0-n)  +  VB(P,^0,n). 
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Fig.  4--Plane  wave  diffraction  by  conducting  wedge. 


If  the  observat-jn  point  is  not  in  the  vicinity  of  a  shadow  boundary 
U*iji0  =  180°)  and  for  p  sufficiently  large,  the  diffracted  field  can 
be  expressed  by  [26] 

h  -jkp 

(19)  H.(p,'>)  £  [D( ,n )  +  D(i^+^0,n)]  — — 

J  p 

where 


cos  J  -  cos  £  | 
n  n  I 

is  the  asymptotic  diffraction  coefficient  in  GTD.  The  total  magnetic 
field  can  then  be  expressed,  using  the  above  form  of  the  diffracted 
field,  by 
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(20) 


H*(p,*)  = 

i  v  p~jkp 

£  +  C(*,«0,n)  5 - 

JT 

where  C  is  independent  of  the  range  (p).  The  total  surface  current 
density  on  a  perfectly  conducting  surface  can  be  obtained  by  applying 
the  boundary  conditions,  and  is  given  by 

J  =  (n  x  z)H^ 

where  n  is  the  unit  outward  normal  vector  to  the  surface.  Using  the 
GTD  solution  for  the  total  magnetic  field,  the  total  surface  current 
density  along  the  wedge  walls  is  given,  for  p  sufficiently  large,  by 

(21)  J6™  =  (n  x  z)  [Hj  +  H[  +  H<] 

*  f  +  Jr  +  (n  x  z)  C(*,*  ,n)  . 

0  JT 

Equation  (21)  is  also  valid  for  wedges  with  local  inhomogenities  around 
the  tip,  assuming  that  the  wedge  surface  is  not  in  the  vicinity  of  a 
shadow  boundary  and  p  is  sufficiently  far  away  from  the  local  in¬ 
homogenities. 

Note  that  C  is  now  assumed  to  be  an  unknown  constant,  although  it 
can  be  obtained  from  a  canonical  solution  for  the  perfectly  conducting 
case.  The  hybrid  approach  can  be  applied  to  a  set  of  problems  where 
the  coefficient  C  is  not  known  at  the  outset. 

The  magnetic  field  integral  equation  for  cylinders  with  arbitrary 
anisotropic  surface  impedance  is  given  by  [29] 


(22)  -hJ(s')  -  }j(s')  +  \~Js(s)  H^)(k|S,-?|)dS 

s  0 

+  ik[j  (s)  .(p H{2>(KD?'  -f|)ds 

Js  |  P ' -p I 

where  the  time  dependence  is  assumed  to  be  e+Jut  and  is  suppressed,  the 
integrals  are  extended  over  the  contour  S  of  the  cylinder  cross  section, 
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n  is  the  outward  normal  to  the  contour  S,  and  ng  is  the  intrinsic  im¬ 
pedance  of  free  space.  Consider  the  impedance-loaded  wedge  structure 
as  illustrated  in  Fig.  5,  the  total  surface  current  density  can  be 


defined  by 


(23)  J  = 


GTD 

h 

h  <  h  <  00 

on 

h-wall 

c  -  - 

MM 

h 

0  <  h  <  h„ 

on 

h-wall 

—  —  c 

GTD 

8 

V 

X 

V 

X 

on 

x-wal 1 

X 

c  -  - 

MM 

X 

0  <  x  <  x„ 

—  —  c 

on 

x-wall 

PTR 

where  J  is  defined  by  Eq.  (21)  which  is  valid  only  away  from  the 
impedance  or  the  tip.  The  method  of  moment  current  (JMM)  around  the 
inhomogenities  can  be  defined  by  simple  basis  functions  such  as 


(24)  JMM  =  I  ampni(s*s  ) 

'  '  ^1  m  m  nr 

where  Pm(s-sm)  are  orthogonal  pulse  functions  with  weight  am. 

Applying  the  point  matching  technique  at  the  midpoint  of  each  of 
the  N  pulse  segment,  one  obtains 


,  ,  / .GTD*  ,  .  /  .GTD, 
Vnn,+  Ln<Jh  >  +  Ln<°x  > 


for  1  <  n  <  N 
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GO 


Fig.  5-- Impedance  loaded  wedge  applied  in  hybrid  methods. 


e  =  2.71828, 


7  =1.781..  is  the  Euler's  constant, 

9n  =  ‘Hz  <Vyn> 


and 
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(26) 


^GTD>*HS 


H<2,<k|v?'|) 


jGTD 

6TD 


+  j 


(?n-?)-n 


h1(2)(Mp„-?'I) 


ds 


Note  that  the  first  term  in  Eq.  (26)  vanishes  because  Zs  =  0  in  the 
6TD  region.  Substituting  Eq.  (21)  into  Eq.  (25)  gives 


(27) 


a  2.  +  L 

m  nm  n 


+  Jh +  “h 


+  L. 


J  +  J  +  D 


1  <  n  <  N 


where  the  second  term  exists  on  the  h-wall  and  the  third  term  on  the 
x-wall.  Since  the  incident  and  reflected  currents  (J1  and  Jr)  can  be 
easily  obtained  by  geometrical  optics,  they  can  be  transferred  to  the 
source  side.  Since  Dv  and  Dh  are  just  conplex  constants,  by  applying 
the  linear  property  of  the  operator  L  ,  Eq.  (27)  can  be  written  as 


(28) 


where 


I  amli 
m 


ru¬ 


in  mn 


DhLn 

n  n 


°xLn 


j 


for  1  <  n  <  N 


=  -Hz<W 


■  Ln<JX>  -  L„<JX> 


(29) 


<k !?„-?!  W 


(pn-tT)’n 


Hl2,(k|Pn-?l) 
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(30) 


e^l  k  fhc 

Ijtt  I  fL 


,(2) 


(k Ip„-p I )+J 


x 


dh' 


H<2)(k|?n-?|)+j 


(31)  Ln(jX)^ 


jkx '  cosiji 

X  e  0  dx'  for  0°  <  *  <  180° 

—  o  — 


0 


otherwise 


h<2)(Mv?IW 


<32>  Ln(JXH 


jkh'cos(  e-180°+ij)  ) 
x  e  0  dh' 


for  180°-e<^o<360°-3 


0 


otherwise  . 


So  far,  we  only  have  N  equations  obtained  by  point  matching  in 
the  moment  method  region.  However,  there  are  N+2  unknowns,  N  coeffi¬ 
cients  of  the  current  pulses  (an's)  in  the  moment  method  region  and 
two  diffraction  coefficients  (Dx  and  Dh),  one  in  each  GTD  region. 


Two  additional  equations  are  obtained  by  introducing  one  addi‘5onal 
matching  point  in  each  of  the  two  GTD  regions.  Matching  at  (x^O)  on 
the  x-wall  gives 
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(33) 


N 

I 

m=l 


a  i  +  D,  L 
in  nm  h  n 


e~jkh 

7T 


j 


+  1  D  - 
2  x 


-jkx^ 


F 


-  -Hz(Xd-°)  -  r  <JX>*=* H-L„K+Jhr> 


Similarly,  matching  on  the  h-wall  gives 


(34) 


N 


l 

m=l 


a  « 
m  mn 


L 


n 


■  -H>=hd>  '  r  <JH>h=hd-L„<JX> 


Now  we  have  N+2  equations  to  solve  for  N+2  unknowns.  This  system 
of  simultaneous  linear  equations  can  be  easily  solved  by  the  numerical 
matrix  inversion  technique. 

The  integrations  over  the  GTD  region  are  infinite  in  theory;  how¬ 
ever,  in  practice,  one  can  approximately  carry  out  the  integration 
over  a  finite  region  because  the  Hankel  function  causes  the  integrand 
to  become  negligible  for  large  value  of  p'. 

The  above  approach  is  valid  except  near  the  shadow  boundaries. 

In  order  to  solve  for  the  diffracted  currents  on  or  near  a  shadow 
boundary  one  should  represent  the  diffracted  currents  in  terms  of  an 
infinite  series  as  given  by  [28] 


J 


d 


e-Jkt 

<j-p)*  • 


However,  in  practice,  if  the  currents  are  analyzed  directly  along  a 
shadow  boundary,  then  the  diffracted  current  can  be  simply  expressed 
by  [28] 
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The  hybrid  method  gives  good  results  for  finite  bodies.  But  it 
fails  to  determine  the  diffracted  field  of  an  infinite  structure 
accurately.  The  values  of  and  Dx  are  found  to  be  sensitive  with 
respect  to  the  position  of  the  GTD  matching  points  in  the  case  of 
infinite  geometries.  This  appraoch  also  reguires  extensive  numerical 
integration  to  evaluate  the  coefficients  L  (JG'D)  and  the  excitation 
g^  of  the  matrix  equation.  n 

C.  A  Modified  Hybrid  Method 

The  key  to  the  modified  hybrid  method  is  to  separate  the  total 
surface  current  of  the  impedance-loaded  wedge  into  two  components, 
namely:  the  current  for  the  perfectly  conducting  case,  Jc:  and  the 
change  in  current  due  to  the  presence  of  the  surface  impedance,  aJ. 
Substituting  J  =  J  +  aJ  into  Eq.  (22)  gives 


(35) 


-h>')-K<s 


)  +f  j  ^JC(S)  H<2’(k|p--p|)ds 

j|^H<2>(k|p’-?l>dS 


^SJ(s)  H<2)(k|?'. 


)ds 


Since  J  satisfies  the  special  case,  Z  =0,  of  the  integral  equation, 
i.e.,  c 

■w  A 

(36)  H^(s')  =  *-Jc(s')  ♦  f  Jc(s)  ^7'“--  H{2)(k|p'-?|)ds 

Js  |p'-p| 

Eq.  (35)  reduces  to 
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(37)  -  | 


j  - - 0  (s)H<2)(ki?'-H)ds  =  t  ‘ J ( s ’ } 

-  s  ’o 


+  I  f  k -  .SJ(s)H<2)(k!?'-?!)<1S  +  iJ(s) 


4  )  s  r'o 


.p  -p)-n 


•  h|2  ^  ( k  I p"'  -p !  )ds  . 


Define  the  change  in  surface  current  density  as 


h  <  h  < 
c  -  - 


on  h-wall 


0  <  h  <  h 


on  n-wall 


•D  =•; 


lsJGTD 


x  <  x  <  <* 
c  -  - 


on  x-wal 


0  <  x  <  x 
-  -  c 


on  x-wal 1 


where  AJ  ^  has  the  form  of  aD  (e  J  h  )//p"»  and  the  moment  method 
incremental  current  can  be  expressed  in  terms  of  simple  orthogonal 
pulse  function,  such  as 


aJ  =  y  :  ,t  P  (s-s  ) 
m  nr  m' 
m=l 


By  point  matching  at  the  center  of  each  pulse,  one  obtains 


N  /  p-jkx 

L  +  A.D  L  I  - - - 

m=l  mm  X"\JT 


AD,  Lj 


h  "n 


h  ; 


=  g  1  <  n  <  N 
n  —  — 


where  a  and  L  are  the  same  as  in  section  B  and 
mn  n 


gn  =  "  A 


's  'o 


Jc(s)H^2)(kj"  -"j)ds  . 
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Again,  two  additional  equations  will  be  obtained  from  the  matching 
in  the  GTD  regions.  A  system  of  N+2  simultaneous  linear  equations 
will  be  solved  for  Act's  and  aD  ,  aD.,  by  numerical  matrix  inversion 
technique.  m  x  n 

It  should  be  noted  that  two  of  the  integrals  in  Eq.  (37)  are 
limited  to  only  the  contour  where  Zs  is  non-vanishing.  Consequently, 
for  conducting  scatterers  with  small  inhonpgenities  with  Zs  t  0 
those  two  integrals  can  be  readily  and  accurately  evaluated.  Further¬ 
more,  as  a  consequence  of  separating  the  current  Jc  for  the  perfectly 
conducting  body,  the  geometrical  optics  currents  are  not  included 
in  the  third  integral  of  Eq.  (37)  which  involves  aJ.  Thus  the  third 
integral  converges  rapidly  because  only  diffracted  components  are 
included  in  t  J.  Those  features  of  the  aJ  approach  result  not  only 
in  significantly  improved  efficiency,  but  also  provide  reasonable 
accuracy  for  the  diffracted  fields.  Computer  solutions  of  the  integral 
equations  typically  result  in  accuracies  of  a  few  percent.  Thus  aJ 
and  the  resulting  diffracted  current  can  be  calculated  to  within  a 
few  percent.  However,  the  original  .ybrid  method  of  Eq.  (22)  yields 
the  same  accuracy  for  the  total  current.  This  results  in  very  poor 
accuracy  for  the  diffracted  component  which  is  typically  a  small 
fraction  of  the  total  current.  This  "aJ"  approach  has  then  a 
distinct  advantage  if  one  wishes  to  isolate  the  diffraction  of  some 
isolated  anomally  in  the  presence  of  an  edge. 
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CHAPTER  IV 

THE  REDUCTION  OF  BACKSCATTERED  FI *LD  —  TWO-DIMENSIONAL  CASE 
A.  Introduction 


Edge  diffracted  fields  become  a  domi na.it  contributor  to  the  scat¬ 
tered  fields  for  a  wide  variety  of  geometries.  These  usually  are 
significant  when  the  observation  point  is  not  illuminated  by  a  specular 
field.  Some  examples  include:  the  cone  for  near  axial  incidence,  the 
leading  or  trailing  edge  of  a  wing  for  sufficiently  low  frequencies. 

These  fields  are  generally  considered  to  be  low  but  this  is  only 
when  considered  with  respect  to  the  specular  contribution  from  a 
surface. 

It  is  the  goal  of  this  study  to  achieve  further  reduction  of 
back  scattered  fields  from  the  edge  diffraction  mechanism  over 
reasonable  range  of  incidence  angles  and  broad  frequency  band  by  use 
of  antenna  impedance  loading.  Absorbing  material  was  often  used  in 
the  study  of  backscatter  reduction.  However,  it  is  usually  necessary 
to  cover  a  substantial  part  of  the  body  of  an  aircraft  or  space  craft. 

The  presence  of  a  large  absorber  could  seriously  affect  the  performance 
of  the  aircraft.  This  is  not  desirable  because  it  would  require  the 
redesign  of  the  whole  aircraft  structure.  In  this  chapter,  the 
possibility  of  designing  an  antenna  to  scatter  energy  at  least  in 
some  direction  or  range  of  directions  in  such  a  way  to  cancel  or 
reduce  the  edge  diffracted  field  in  that  range  of  direction  is  con¬ 
sidered.  Surface  impedance  is  used  to  model  slot  antenna.  The  hybrid 
methods  discussed  in  Chapter  III  are  employed  to  treat  various  geometries 
and  the  results  for  optimum  loading  are  also  presented. 

B.  The  Triangular  Cylinder 

The  reaction  integral  equation  and  the  Galerkin's  method  with  a 
piece-wise  sinusoidal  basis,  developed  by  Richmond  [30],  is  used  to 
study  the  scattering  by  a  two-dimensional,  triangular  wedge-shaped 
cylinder.  The  scattering  geometry  of  a  triangular  cylinder  is  quite 
suitable  for  providing  information  on  the  effectiveness  of  impedance 
loading  for  backscatter  reduction.  Furthermore,  this  geometry  provides 
efficient  numerical  results.  Results  on  optimum  loading  configurations 
obtained  from  this  geometry  will  be  directly  applicable  to  a  wide 
variety  of  scattering  geometries  such  as  conical  shapes  and  ogive 
shape  or  similar  wing  structures. 

In  order  to  isolate  the  effect  of  the  impedance  loading  a  source 
excitation  is  being  used  that  consists  of  two  line  sources  located  on 
the  cylinder  surface  and  phased  to  null  out  the  effect  of  the  other 
edge.  The  wedge-shaped  cylinder  without  impedance  loading  can  readily 
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be  analyzed  by  the  GTD  or  wedge  diffraction  approach.  The  wedge  dif¬ 
fraction  analysis  for  the  unloaded  cylinder  was  programmed  as  a  check 
case.  The  value  of  the  wedge  diffraction  check  was  demonstrated  when 
it  revealed  an  insufficient  number  of  piece-wise  sinusoidal  monopoles 
were  used  in  the  reaction  integral  equation  program.  The  results  ob¬ 
tained  by  using  4  segments  in  one  wavelength  and  the  GTD  results  are 
shown  in  Fig.  6.  The  difference  between  these  two  results  is  sig¬ 
nificant.  However,  if  8  segments  in  one  wavelength  are  used,  the 
result  agrees  well  with  the  GTD  result  (see  Fig.  7). 

Due  to  limitations  on  the  dimensions  of  the  wedge  cylinder  that  can 
practically  be  used  in  the  reaction  integral  equation  program  on  the 
DATACRAFT  6024  computer,  results  can  only  be  obtained  for  cylinders 
on  the  order  of  2.5A  or  smaller  in  width.  However,  several  promising 
results  have  been  obtained  using  the  integral  equation  program.  These 
results  indicate  that  considerable  reduction  can  be  achieved  in 
the  backscattering  from  the  edge  diffraction  mechanism.  In  one  case 
a  reduction  of  14  dB  was  obtained  using  a  real  impedance  loading  of 
377  ohms  located  right  at  the  edge  as  shown  in  Fig.  8.  The  com¬ 
puted  results  also  show  that  similar  reduction  can  be  achieved  by  a 
reactive  loading  of  Z  =  -jl 88  ohms  located  5/8x  from  the  edge,  as 
shown  in  Fig.  9.  s 

The  backscattered  fields  of  such  a  wedge  geometry  have  also  been 
computed  as  a  function  of  incidence  angle  and  surface  impedance.  This 
analysis  shows  the  reduction  of  the  edge  diffracted  field  discussed 
above  is  maintained  from  angles  of  zero  (preceding  case)  to  about  60° 
when  the  surface  impedance  is  at  the  edge.  Figure  10  shows  the  back- 
scattered  fields  with  zero  and  188  ohms  surface  impedance.  Note  that 
the  lobes  are  decreased  from  40°-100°.  For  plane  wave  incidence  as 
shown  in  Fig.  10  the  effect  of  the  edge  that  does  not  have  impedance 
loading  cannot  be  isolated  from  the  edge  with  loading.  Consequently, 
the  reduction  in  backscatter  as  shown  in  Fig.  10  represents  a  considerable 
reduction  achieved  by  the  loaded  edge.  This  corresponds  to  the  re¬ 
duction  achieved  in  Fig.  8,  but  also  demonstrates  the  wide  angle 
coverage  that  can  be  achieved.  Thus,  the  above  analysis  does  indicate 
that  the  desired  control  can  be  achieved. 

C.  The  Infinite  Wedge 

The  hybrid  methods  discussed  in  Chapter  II  were  programmed  for 
the  analysis  of  the  infinite  wedge  with  impedance  loading.  The  total 
current  on  a  perfectly  conducting  wedge  surface  was  calculated  by  the 
hybrid  method  of  Eq.  (22).  The  current  was  sampled  using  a  0.1X  o** 
smaller  spacing  between  sampling  points  in  the  MM  region.  The  moment 
method  region  is  lx  in  length  and  each  GTD  matching  point  is  1.25  a 
away  from  the  edge  tip.  Results  for  the  surface  current  density  and 
the  diffracted  field  of  a  30°  conducting  wedge  with  plane  wave  in¬ 
cidence  are  compared  with  GTD,  as  shown  in  Figs.  11-12.  Although  the 
original  hybrid  method  was  found  to  yield  reasonable  results  for  the 
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Fig.  6--Comparison  of  patterns  from  GTD  and  integral 
equation  solution. 


24 


Fig.  8—Comparison  of  radiation  patterns  with  and  without 
surface  impedance. 


26 


Fig.  9--Compari son  of  radiation  patterns  with  and  without 
surface  impedance. 
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Fig.  10— Plane  wave  backscatter  from  loaded  and  unloaded 
triangular  cylinder. 
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total  current  and  high-level  fields  of  the  scatterer,  it  fails  to 
accurately  predict  low-level  fields  such  as  diffracted  fields,  as 
seen  in  Figs.  11-12.  Consequently,  the  original  hybrid  method  was  found 
to  be  inadequate  to  analyze  the  performance  of  impedance  techniques. 
After  some  research  it  was  found  that  the  original  hybrid  method 
could  be  extended  to  provide  accurate  results  for  low  level  dif¬ 
fracted  fields. 

As  has  been  noted  the  key  to  the  modified  hybrid  method  is  to 
seperate  the  surface  current  of  an  impedance  loaded  scatterer  into 
two  components:  the  current  Jc  for  the  perfectly  conducting  scatterer 
and  the  change  in  current  aJ  due  to  impedance  loading.  Since  the  cur¬ 
rent  Jc  for  the  perfectly  conducting  wedge  is  known,  the  current  aJ 
on  an  impedance  loaded  wedge  can  be  solved  using  the  modified  hybrid 
method.  The  current  aJ  is  sampled  using  a  O.U  or  smaller  spacing 
between  sampling  points  in  MM  region.  Since  the  current  aJ  decays 
rapidly  away  from  any  impedance  loadings,  the  integrals  of  the  modi¬ 
fied  method  converge  rapidly.  This  results  in  improved  efficiency  as 
well  as  accuracy  compared  to  that  of  the  original  hybrid  method. 

The  results  of  reduction  in  backscatter  that  can  be  achieved  by 
impedance  loading  are  shown  in  Figs.  13-18.  A  reduction  of  more  than 
15  dB  is  obtained  in  an  angular  range  of  20°  from  grazing  incidence  for 
both  30°  and  60°  wedge.  These  reductions  were  also  demonstrated  at 
three  different  frequencies  representative  of  a  3:1  frequency  band. 

The  impedance  loading  required  over  the  continuous  frequency  range 
was  also  calculated.  The  required  surface  impedance  loading  over 
the  3:1  frequency  range  and  the  reduction  in  backscatter  from  the  30° 
infinite  wedge  are  given  in  Table  I .  The  surface  impedance  is  plotted 
on  a  Smith  chart  in  Fig.  19. 

D.  The  Square  Cylinder 

The  accuracy  of  the  modified  hybrid  method  was  checked  against  the 
integral  equation  or  moment  method  solution  for  the  4a  square  cylinder 
with  impedance  loading  as  shown  in  Fig.  20.  This  cylinder  size  is 
close  to  the  limit  possible  with  the  pure  moment  method.  The  accuracy 
of  the  modified  hybrid  method  was  found  to  be  better  than  6%  even  for  low 
level  fields.  The  computed  results  of  the  two  methods  are  so  close  that 
the  difference  can  not  be  seen  on  Fig.  20. 

The  scattered  field  from  a  square  cylinder  is  shown  in  Fig.  21. 

The  results  for  this  case  were  calculated  by  all  three  methods:  the 
moment  method,  the  modified  hybrid  method  and  the  original  hybrid 
method.  It  was  found  that  the  original  hybrid  method  gives  good 
results  for  the  4x  cylinder  in  contrast  to  the  case  of  the  infinite 
wedge.  Thus  it  can  be  concluded  that  the  original  hybrid  method  is 
sufficiently  accurate  for  finite  size  scatterers  for  which  the  in¬ 
tegrations  involved  in  the  original  hybrid  method  are  not  truncated. 
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Fig.  1 8 — Edge  diffracted 


TABLE  I 

REDUCTION  OF  BACKSCATTER  AT  GRAZING  INCIDENCE  FOR  30°  WEDGE 


Vno 

d/x 

w/a 

Reduction 

1.75+jl.6 

0.05 

0.05 

25 

1.85+jl.O 

0.06 

0.06 

23 

1.7  +j0.3 

0.07 

0.07 

19.5 

1. 45- jO. 1 

0.08 

0.08 

18 

1.15-jO. 3 

0.09 

0.09 

21 

0.9  -j0.5 

0.1 

0.1 

30 

0.7  -jO.55 

0.11 

0.11 

21.5 

0.5  -j0.6 

0.12 

0.12 

21 

0 . 35- j  0.6 

0.13 

0.13 

23 

0.23-jO. 6 

0.14 

0.14 

27 

0.1  -j0.6 

0.15 

0.15 

24 
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Fig.  19--Surface  impedance  to  reduce  backscatter 
from  an  infinite  wedge  over  a  3:1 
frequency  band. 
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Fig.  21— Scattered  field  of  conducting  square  cylinder 
witii  surface  impedance  loading. 
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E.  A  Two-dimensional  Wing  Model 


The  original  hybrid  method  was  used  to  analyze  the  two-dimensional 
wing  model  shown  in  Fig.  22(a).  This  wing  model  should  provide  a  good 
simulation  for  a  practical  aircraft  wing  cross  section.  However,  the 
small  wedge  (angle  %6°)  of  the  wing  model  of  Fig.  22(a)  was  found  to 
cause  a  computational  problem  in  the  application  of  the  hybrid  method. 
A^ter  investigation  of  this  problem  it  was  discovered  that  the  small 
weuge  angle  causes  some  of  the  matching  points  used  in  the  point 
matching  region  around  the  edge  to  be  too  close  for  practical  computation 
as  shown  in  Fig.  22(b).  This  problem  was  overcome  by  truncation  of  the 
wedge  tip  as  shown  in  Fig.  22(c).  The  truncation  still  permits  a 
practical  analysis  of  the  edge  loading  because  the  trailing  edge  of 
the  wing  is  very  thin  in  wavelengths. 

The  current  in  the  MM  regions  is  again  sampled  using  a  O.U  or 
smaller  spacing  between  sampling  points.  Since  the  primary  interest 
is  in  the  angular  range  which  is  close  to  grazing  to  the  wedge  arm, 
the  diffracted  currents  from  the  wedge-cylinder  junction  are  analyzed 
along  or  near  the  shadow  boundary.  The  ^orm  D0e“Jk£  +  D]  (e_J^^)//T 
is  used  for  the  diffracted  currents  from  the  wedge-cylinder  junction, 
while  D  (e*Jkh)//F  is  used  for  the  diffracted  currents  from  the  wedge 
tip,  as  shown  in  Fig.  23. 

The  backscatter  reduction  of  the  two-dimensional  wing  model  as  a 
function  of  wing  size  and  geometry,  including  the  effect  of  the  leading 
edge,  has  been  investigated.  The  scattered  fields  for  Zs  =  0  and 
Zs  =  n0  on  the  leading  edge  were  checked  by  using  the  hybrid  method 
and  the  moment  method  for  a  -  0.075a,  0.1a,  0.15a,  0.2a,  0.25a  and 
0.3A.  Figures  24-25  show  some  of  the  typical  checks.  Tables  II  and 
III  show  the  backscatter  reduction  from  the  leading  and  trailing 
edge  for  various  wing  size.  The  leading  edge  is  dominant  for  wing 
radii  greater  than  0.1  A.  However,  the  leading  edge  contribution  can 
be  reduced  by  more  than  12  dB  for  a  ^0.2A,  as  shown  in  Fig.  26.  This 
is  achieved  by  covering  the  leading  edge  with  a  surface  impedance 
Zs  =  no-  Eor  wing  radii  smaller  than  O.lA,  the  trailing  edge  con¬ 
tribution  is  more  important.  By  loading  the  trailing  edge  with  a  0.1a 
wide  surface  impedance,  the  backscatter  for  a  thin  wing  can  be  reduced 
by  10.19  dB,  as  shown  in  Fig.  27.  The  backscatter  reduction  from  the 
leading  and  trailing  edge,  respectively,  as  a  function  of  incidence 
angle  is  also  calculated  and  shown  in  Tables  IV  and  V.  It  is  seen 
from  Fig.  28  that  further  reduction  can  be  obtained  if  both  edges  are 
loaded.  The  reduction  in  Fig.  28  is  11.33  dB,  compared  to  5.79  dB 
when  only  the  leading  edge  is  loaded  and  5.37  dB  when  only  the  trailing 
edge  is  loaded. 
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Fig.  22— Point-matching  problem  for  small-angle  wedge. 
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Fig.  23— MMGTD  method  applied  to  the  two-dimensional 
wing  model . 
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Fig.  24 --Comparison  of  bistatic  scattered  fields  obtained 
using  MM  and  MMGTO  solutions  for  a  wing  foil  with 
impedance  loading. 
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Fig.  25 --Comparison  of  bistatic  scattered  feilds  obtained  using  MM 
and  MMGTD  solutions  for  a  wing  foil  with  impedance 
loading. 
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TABLE  II 


BACKSCATTER  REDUCTION  BY  LOADING  THE  LEADING  EDGE 

Radius  (a)  Backscatter  Reduction  (dB) 

Moment  Method  MMGTD 


0.3A 

17.77 

18.00 

0.25A 

14.89 

14.92 

0.2A 

12.28 

12.24 

0.1 5A 

9.42 

9.40 

0 . 1  A 

5.85 

5.79 

0.075A 

3.30 

3.30 

TABLE  III 

BACKSCATTER  REDUCTION  BY  LOADING  THE  TRAILING  EDGE 

Radius  (a)  Backscatter  Reduction  (dB) 

0.3* 

0.25a 
0.20A 
0.15a 
0.10a 
0.075a 


d  *  w  =  0.  IX 


IT  *  1-0- j-0.6 
#0 


-1.0 

0.05 

1.29 

2.06 

5.37 

10.19 
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Fig.  26--Comparison  of  bistatic  scattered  fields  for  wings 
with  and  without  impedance  loading. 
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Fig.  27— Comparison  of  bistatic  scattered  fields  for  wings  with 
and  without  impedance  loading. 
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TABLE  IV 


BACKSCATTER  REDUCTION  OF  WING  MODEL  AS  A  FUNCTION  OF  INCIDENCE  ANGLE 


Incidence  Angle  (e) 


Backscatter  Reduction  (dB) 


0 

3 

8 

10 


18.00 

17.57 

16.47 

12.08 


TABLE  V 

BACKSCATTER  REDUCTION  OF  WING  MODEL  AS  A  FUNCTION  ON  INCIDENCE  ANGLE 


Incidence  Angle  (0) 


Backscatter  Reduction  (dB) 


0 

5 

10 


10.19 

11.13 

3.12 


d  *  w  =  0.1  X 
^7  *  I.O-j-O.6 
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Fig.  28-Comparison  of  bistatic  scattered  fields  for  wings 
with  and  without  impedance  loading. 
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F. 


The  Thin  Strip 


The  infinitely  long  thin  strips  of  1  and  2  wavelength  wide  and  0.1 
wavelength  thick  are  studied  by  moment  method.  The  backscatter  reduction 
will  be  used  to  check  the  results  obtained  from  a  three-dimensional 
finite  thin  plate  discussed  in  Chapter  V.  Typical  results  for  1  and 
2  wavelength  strips  loaded  with  a  O.lx  wide  surface  impedance  are 
shown  in  Figs.  29  and  30.  Tables  VI  and  VII  show  the  backscatter 
reduction  as  a  function  of  incidence  angle. 
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Fig.  29--Comparison  of  oistatic  scattered  fields  for  strips 
with  and  without  impedance  loading. 
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Fig.  30- -Comparison  of  bi static  scattered  fields  for  strips 
with  and  without  impedance  loading. 
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TABLE  VI 


BACKSCATTER  REDUCTION  OF  THIN  STRIP  AS  A  FUNCTION  OF  INCIDENCE  ANGLE 


Incidence  Angle  (e) 


Backscatter  Reduction  (dB) 


0 

3 

5 

10 


15.61 

14.97 

13.83 

8.97 


TABLE  VII 

BACKSCATTER  REDUCTION  OF  THIN  STRIP  AS  A  FUNCTION  OF  INCIDENCE  ANGLE 


Incidence  Angle  (e)  Backscatter  Reduction  (dB) 


0 

3 

5 

10 


16.45 

15.89 

13.39 

5.45 


0.1  X 

d*w«O.IX  —  «  0.7 -j- 0.6 

Vo 
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CHAPTER  V 

THE  REDUCTION  OF  BACKSCATTERED 
FIELD  —  THREE-DIMENSIONAL  CASE 


A.  Introduction 

The  results  previously  discussed  yield  the  necessary  surface 
impedance  loadings  required  to  achieve  reduction  in  the  backscatter. 
However,  a  way  to  practically  implement  the  required  surface  impedance 
is  also  needed.  One  of  the  most  promising  ways  to  implement  the  sur¬ 
face  impedance  is  a  col  linear  array  of  rectangular  slots  located  along 
the  edge. 

In  this  chapter,  the  thin  conducting  plate  with  a  loaded  slot  is 
treated.  The  reductions  in  backscatter  of  various  sizes  of  plates  are 
also  illustrated.  It  is  concluded  that  the  backscatter  reduction  can 
be  achieved  over  a  2:1  frequency  band  in  a  wide  angular  aspect. 

B.  The  Thin  Plate 


Since  the  thin  square  plate  with  a  loaded  slot  antenna  can  provide 
the  information  about  the  affect  of  loading,  the  backscatter  reduction 
of  the  plate  is  investigated.  For  near  grazing  angles  of  incidence, 
the  trailing  edge  of  the  plate  is  the  dominant  contributor  to  the 
scattering  of  the  component  of  the  electric  field  which  is  normal  to 
the  edge  (vertical  polarization).  On  the  other  hand,  the  leading  edge 
contribution  is  the  dominant  factor  to  the  scattering  of  the  component 
of  the  electric  field  which  is  parallel  to  the  edge  (parallel 
polarization). 

The  reduction  of  the  backscatter  of  the  plate  for  the  vertical 
polarization  can  be  achieved  by  inserting  a  slot  sufficiently  near  the 
trailing  edge  of  the  plate.  The  axis  of  the  slot  should  be  oriented 
parallel  to  the  trailing  edge  in  order  for  the  slot  to  have  maximum 
effect.  A  lumped  load  impedance  may  be  inserted  in  the  slot  to  tune 
the  slot.  The  size,  location  and  the  load  impedance  of  the  slot  have 
to  be  chosen  properly  in  order  to  achieve  a  reduction  in  backscattering 
cross  section  over  a  wide  frequency  band. 

The  backscattering  analysis  for  the  plate  with  loaded  slot  has 
been  carried  out  by  approximating  the  conducting  plate  with  a  slot  by 
a  conducting  wire  grid  model.  The  integral  equations  for  the  scat¬ 
tering  of  electromagnetic  plane  wave  by  the  wire  grid  structure  are 
then  solved  by  the  numerical  techniques.  A  typical  wire  grid  model  of 
a  half-wave  slot  in  one  wavelength  square  plate  is  illustrated  in  Fig. 
31.  The  lumped  load  impedance  of  the  slot  antenna  is  obtained  by 
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inserting  a  proper  terminal  load  impedance  at  the  center  c  of  the  wire 
segment  abc  in  Fig.  31 . 


r 

r 

0 

« 

>c 

b 

Fig.  31 --Conducting  wire  grid  model  of  a  half-wavelength 
long  slot  in  a  one  wavelength  square  con¬ 
ducting  plate. 


The  wire-grid  structure  is  considered  as  an  array  of  thin  wire 
segments.  Current-sampling  points  or  terminals  are  defined  at  each 
wire  junction. 

Assume  that  each  segment  in  the  wire  structure  has  a  circular 
cylindrical  surface  and  the  wire  radius  "a"  is  much  smaller  than  the 
wavelength  a,  and  the  wire  length  is  much  greater  than  the  radius. 
Furthermore,  the  circumferential  component  Ja  of  the  surface  current 
density  will  be  neglected.  If  we  consider  tne  axial  component  of  the 
surface  current  density  to  be  independent  of  <f>,  then  the  surface  cur¬ 
rent  density  of  the  thin  wire  structure  is  given  by 


(40) 


iUBl  -  till 

2-rra  2tt3 


where  £  is  the  coordinate  along  the  wire  axis  and  I ( it)  is  the  total 
current.  Substituting  Eq.  (40)  into  Eq.  (8)  and  using  the  impedance 
boundary  condition  of  Eq.  (9)  leads  to 
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(41)  IC)  (Ej-Z^JdH  =  Vm 

where  L  denotes  the  overall  wire  length  and 


(42)  Vm  =  fff(Ji-Em-Mi.HTn)dv 

J 


(43) 


.  1_  f2lT 

'  2-  Jn 


i-f' 


(44) 


um  -  1 

%  -  27  j 


2tt 


0 


i-  H™  d$ 


In  order  to  solve  for  the  unknown  current  distribution  in  Eq.  (41), 
a  sinusoidal  dipole  is  used  as  a  test  source.  The  current  distribution 
on  a  V  test  dipole  is 


(45) 


£-|P.  sinh  y(£-£.j  ) 
sinh  y«~| 


P2  sinh  y(£2-£) 
sinh  y&2 


where  P1  and  P2  are  unit  pulse  functions  and 

Y  =  ji ojue 


The  current  distribution  on  the  wire  structure  is  expanded  in  a  finite 
series  of  sinusoidal  functions, 


(46) 


1(4) 


N 

"A 


tnU) 


Substituting  Eq.  (46)  into  Eq.  (41)  gives 


(47) 


N 

y  i  z  =  v 

n  mn  m 


m  =  1,2  N 
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where 


(48) 


F„(*>  * 


Equation  (47)  can  be  solved  for  the  unknown  current  I  by  numerical 
techniques.  n 

If  a  lumped  impedance  Z,n  is  inserted  at  the  m*h  current  sampling 
points  of  the  wire  structure,  there  is  a  voltage  drop  ImZm  across  the 
load.  Therefore,  in  the  delta  gap  model,  the  lumped  loads  simply 
introduce  a  new  term  in  Eq.  (47)  so  the  right-hand  side  becomes 
Vm~ImZm’  Transfemn9  the  last  term,  we  obtain 

(49)  j,  ‘n(Z™  *  V.J  '  Vn  ' 

Once  the  current  distribution  is  solved,  the  scattered  field  can  be 
readily  calculated. 

The  main  computer  program  for  the  analysis  of  the  backscattering 
of  a  plate  is  based  on  the  computer  subroutines  developed  by 
Richmond  [33].  In  order  to  model  the  plate  with  slot  accurately, 
it  is  found  that  the  spacing  between  wire  grid  should  be  as  small  as 
a/8,  at  least  in  the  vicinity  of  the  slot. 

The  plate  with  a  slot  can  be  considered  as  an  antenna.  Using  the 
compensation  theorem  and  the  superposition  principle,  the  scattered 
field  from  a  loaded  antenna  may  be  interpreted  as  a  load  independent 
term  called  the  structural  scattering  plus  a  load  dependent  term 
called  the  antenna  mode  scattering  (see  Appendix  A). 

The  backscattering  patterns  in  the  plane  perpendicular  to  the 
slot  axis  are  calculated  for  different  sizes  of  plates  and  slots  using 
various  load  impedances.  From  Eq.  (52)  (see  Appendix  A)  we  see  that 
if  the  short  and  open  circuit  scattered  field  from  a  loaded  antenna 
are  known  (either  measured  or  calculated)  the  scattered  field  for 
arbitrary  load  can  be  readily  calculated.  This  is  brought  to  the 
author's  attention  by  Richmond  [401.  Considerable  amount  of  computer 
time  can  be  saved  by  using  Eq.  (52)  to  find  the  optimum  load  impedance 
for  backscatter  reduction.  The  backscattered  field  calculated  by 
directly  using  the  wire-grid  model  is  compared  with  that  obtained  by 
Eq,  (52),  as  shown  in  Table  VIII  and  IX.  The  reductions  of  back¬ 
scatter  are  shown  in  Figs.  32-37.  The  load  impedance  which  is  required 
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to  reduce  the  backscatter  from  the  trailing  edge  of  the  square  plate 
over  a  2:1  frequency  band  is  plotted  on  a  Smith  chart  in  Fig.  38.  It 
is  seen  from  Figs.  32-37  that  a  substantial  reduction  (>  15  dB)  in  the 
backscattering  cross  section  can  be  achieved  when  the  slot  is  properly 
loaded.  In  order  to  achieve  backscattering  reduction  for  near  grazing 
incidence  on  either  edge,  two  slots,  one  on  each  edge,  were  used  on  the 
square  plate.  It  is  seen  from  Fig.  37  that  the  reduction  obtained  by 
using  two  slots  is  nearly  as  good  as  that  for  one  slot. 
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TABLE  VIII 

COMPARISON  OF  THE  SCATTERED  FIELD  OF  A  LOADED  SQUARE  PLATE 


0 . 5 x  x  0 . 5x  Plate,  \  =  14.87  -  j - 200 


Calculated  by 
Wire-Grid  Model 

Real  Part  Imaginary  Part 


Calculated 
Eq.  (52) 


0 

0.0 

0.0 

0 

5 

-0.00°99 

-0.00165 

” 

10 

-0.00413 

-0.00652 

- 

15 

-0.00994 

-0.01428 

- 

20 

-0.01924 

-0.02430 

25 

-0.03300 

-0.03553 

30 

-0.05215 

-0.04635 

- 

35 

-0.07730 

-0.05449 

40 

-0.10832 

-0.05709 

45 

-0.14388 

-0.05094 

“ 

50 

-0.1811  o 

-0.03302 

“ 

55 

-0.21556 

-0.00130 

60 

-0.24136 

0.04435 

“ 

65 

-0.25229 

0.10127 

/0 

-0.24316 

0.16376 

75 

-0.21145 

0.22357 

80 

-0.15863 

0.27131 

“ 

85 

-0.09059 

0.29873 

— 

90 

-0.01672 

0.30089 

95 

0.05215 

0.27770 

100 

0.10646 

0.23405 

105 

0.14020 

0.17843 

no 

0.15206 

0.12055 

115 

0.14515 

0.06883 

120 

0.12536 

0.02868 

125 

0.09940 

0.00193 

130 

0.07303 

-0.01260 

135 

0.05005 

-0.01787 

140 

0.03220 

-0.01732 

145 

0.01958 

-0.01 395 

150 

0.01133 

-0.00987 

155 

0.00627 

-0.00626 

160 

0.00329 

-0.00354 

165 

0.00159 

-0.00175 

170 

0.00064 

-0.00069 

175 

0.00015 

-0.00016 

Real  Part 


.00098 

.00412 

.00992 

.01925 

.03303 

.05212 

.07731 

.10833 

.14390 

.18110 

.21549 

.24135 

.25226 

.24306 

.21134 

.15851 

.09049 

.01663 

.05221 

.10649 

.14009 

.15193 

.14500 

.12519 

.09920 

.07283 

.04988 

.03208 

.01954 

.01128 

.00620 

.00328 

.00160 

.00062 

.00012 


by 


Imaginary  Part 

0.00000 

-  .00168 

-  .00654 

-  .01428 

-  .02431 

-  .03555 

-  .04638 

-  .05453 

-  .05711 

-  .05091 

-  .03298 

-  .00128 
.04437 
.10131 
.16377 
.22359 
.27129 
.29861 
.30077 
.27751 
.23386 
.  1  7824 
.12033 
.06867 
.02857 
.00192 

-  .01260 

-  .01781 

-  .01725 

-  .01387 

-  .00977 

-  .00623 

-  .00352 

-  .00171 

-  .00067 

-  .00019 


61 


TABLE  IX 


COMPARISON  OF  THE  SCATTERED  FIELD  OF  A  LOADED  SQUARE  PLATE 


0.6A  x  0.6A  Plate,  =  56.16  -  j . 300 


Calculated  by 
Wire-Grid  Model 


0  Real  Part  Imaginary  Part 


Calculated  by 
Eq.  (52) 


Real  Part 


Imaginary  Part 


0 

0.0 

0.0 

0.00000 

0.00000 

5 

-0.00117 

-0.00207 

-  .00117 

-  .00206 

10 

-0.00489 

-0.00815 

-  .00488 

-  .00817 

15 

-0.01180 

-0.01790 

-  .01179 

-  .01788 

20 

-0.02291 

-0.03059 

-  .02292 

-  .0305  8 

25 

-0.03948 

-0.04503 

-  .03948 

•  .04502 

30 

-0.06282 

-0.05934 

-  .06282 

-  .05936 

35 

-0.09400 

-0.07082 

-  .09400 

-  .07082 

40 

-0.13337 

-0.07582 

-  .13336 

-  .07583 

45 

-0.17991 

-0.06990 

-  .17992 

-  .06990 

50 

-0.23051 

-0.04823 

-  .230  5  2 

-  .04822 

55 

-0.27933 

-0.00665 

-  .27934 

-  .00663 

60 

-0.31763 

0.05672 

-  .31762 

.05675 

65 

-0.33470 

0.13953 

-  .334  6  5 

.13950 

70 

-0.32014 

0.23350 

-  .32011 

.23345 

75 

-0.26748 

0.32433 

-  .26747 

.32423 

80 

-0.17801 

0.39394 

-  .17806 

.39378 

85 

-0.06290 

0.42539 

-  .06305 

.42521 

90 

0.05809 

0.40886 

.05783 

.40870 

95 

0.16177 

0.34612 

.16140 

.34607 

100 

0.22891 

0.25106 

.22845 

.25113 

105 

0.25057 

0.14527 

.2501  3 

.14555 

110 

0.23041 

0.05077 

.23007 

.05121 

115 

0.18207 

-0.01702 

.18190 

-  .01646 

120 

0.12317 

-0.05287 

.12320 

-  .05228 

125 

0.06903 

-0.06094 

.06920 

-  .06040 

130 

0.02873 

-0.05091 

.02902 

-  .05048 

135 

0.00458 

-0.03327 

.00492 

-  .03297 

140 

-0.00609 

-0.01603 

-  .00575 

-  .01586 

145 

-0.00809 

-0.00345 

-  .00780 

-  .00337 

150 

-0.00601 

0.00351 

-  .00578 

.00351 

155 

-0.00310 

0.00587 

-  .00293 

.00583 

160 

-0.00094 

0.00536 

-  .00083 

.00531 

165 

0.00008 

0.00362 

.00013 

.00360 

170 

0.00026 

0.00177 

.00029 

.00176 

175 

0.00010 

0.00046 

.00009 

.00045 

62 


0  25  50  75  100  125  150  175 


9  (degrees) 


fig.  33--Use  of  a  loaded  slot  for  reduction  of  an  edge  diffracted 
field  from  a  flat  plate. 


Fig.  34--Use  of  a  loaded  slot  for  reduction  of  an  edge  diffracted 
field  from  a  flat  plate. 
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Fig.  35--Use  of  a  loaded  slot  for  reduction  of  an  edge  diffracted 
field  from  a  flat  plate. 
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Fig.  36— Use  of  a  loaded  slot  '<  r  reduction  of  an  edge  diffracted 
field  from  a  flat  plate. 
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Fig.  37.  Backscattering  cross  section  (o0e/A  for  the  0  polarization) 
of  a  square  conduction  plate  with  one  and  two  slots. 
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I.OX 


Fig.  38--Load  impedance  required  for  reduction  of  wedge 
diffracted  backscattered  fields. 


CHAPTER  VI 
CONCLUSIONS 


The  application  of  impedance  loading  techniques  for  the  reduction 
of  backscatter  from  radar  targets  were  studied.  The  hybrid  methods 
which  combine  the  GTD  and  the  moment  method  were  used  to  analyze  the 
backscattering  from  various  geometries  such  as  the  impedance-loaded 
infinite  wedge,  thin  strip  and  a  two-dimensional  wing  model.  The 
surface  impedance  was  used  to  simulate  a  slot  antenna  loaded  with  a 
terminal  impedance.  The  original  hybrid  method  was  found  inadequate 
to  analyze  the  reduction  in  scattering  by  edge  loading  an  infinite 
wedge.  A  solution  to  this  difficulty  was  attained  by  the  development 
of  the  aJ  method. 

The  effect  of  the  edge  diffraction  sources  can  be  considerably 
reduced  over  a  3:1  frequency  band  by  the  use  of  a  loaded  slot  near 
or  at  an  edge.  These  reductions  are  found  to  cover  an  angular  range 
of  20°  from  grazing. 

In  order  to  include  the  effect  of  the  leading  edge,  a  two- 
dimensional  wing  model  which  provides  a  good  simulation  of  a  practical 
aircraft  wing  cross  section  was  also  studied.  The  backscatter  re¬ 
ductions  for  the  wing  model  were  calculated  as  a  function  of  wing  size 
and  geometry. 

The  slot-loaded  square  plate  which  is  a  basis  of  an  array  of 
finite  slots  along  the  edge  of  a  three-dimensional  wing  model  was 
analyzed  by  a  wire  grid  model.  The  array  of  loaded  slots  is  one  of 
the  most  promising  ways  to  implement  the  necessary  surface  impedance 
loadings  required  to  achieve  reductions  in  the  backscatter.  The  back¬ 
scatter  reduction  of  the  square  plate  can  be  achieved  at  least  over 
a  2:1  frequency  band.  It  is  also  demonstrated  that  the  use  of  two 
slots  on  the  plate  can  reduce  the  backscattering  for  near  grazing 
incidence  on  either  edge. 
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APPENDIX  A 

THE  SCATTERED  FIELD  OF  AN  ANTENNA 
AS  A  FUNCTION  OF  LOAD  IMPEDANCE 


The  scattered  field  from  an  impedance-loaded  object  can  be 
expressed  as  a  superposition  of  the  field  scattered  by  the  un¬ 
loaded  object  and  the  field  reradiated  by  the  load.  In  attempting 
to  predict  the  scattered  field  from  an  antenna  as  a  function  of  the 
load  impedance  many  researchers  have  prtsented  different  functional 
forms  [8,35,36,37,38,39].  In  the  form  discussed  here,  the  scattered 
field  can  be  interpreted  as  a  load  independent  term  called  the 
structural  scattering  plus  a  term  called  antenna  mode  scattering. 

Consider  the  scattering  of  electromagnetic  wave  by  an  impedance- 
loaded  antenna  as  shown  in  Fig.  39.  Using  the  compensation  theorem, 
the  load  may  be  replaced  by  a  voltage  source.  By  superposition 
principle,  the  scattered  field  is  broken  into  two  parts,  one  due  to 
the  incident  field  on  the  shorted  antenna  and  the  other  to  the 
equivalent  voltage  source  exciting  the  antenna.  Therefore,  the 
scattered  field  by  an  antenna  as  a  function  of  the  load  impedance 
is  given  by 


(50)  E(Z  )  =  E(0)  -  I(0)Er 

£  a 

where 

E (0 )  =  the  field  scattered  by  the  antenna  with  Z  =0,  that  is, 
a  short  circuit  for  the  load, 

1(0)  =  the  antenna  terminal  current  with  =  0, 

E  =  the  field  radiated  by  the  antenna  when  excited  by  a 
unit  current  source, 

and 

Z, ,Z  are  antenna  and  load  irpedance,  respectively. 

a  Xr 

Let  the  load  impedance  approach  infinity,  that  is,  open  circuit, 
then  Eq.  (50)  becomes 

E(»)  =  E(0)  -  I(0)F 

or 

(51)  I(0)F  =  E(0)  -  #(»)  . 


71 


E1 


I(0)ZjZa 


Z*+Zfl 


+ 


IU0)2j 

h+Z0 


Fig.  39— Compensation  and  superposition  theorem  applied 
to  scattering  problem. 


Substituting  Eq.  (51)  into  Eq.  (50)  gives 

Z£[E(0)  -  E(-)] 


(52) 


E(Zt)  =  E(0)  - 


VZa 


Za  E(0)  +  7l  E(») 

S'^rrra 


Equation  (52)  can  be  used  to  calculate  the  scattered  field  for  arbitrary 
load  if  the  open  and  short  circuit  scattered  field  from  the  antenna  are 
known. 
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